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Abstract 

We propose a hybrid rates-credit model for pricing credit default swaps (CDS) in the presence of a 
correlation between rates and credit. We suppose a beta-blend model for rates and a Black-Karasinski 
model for credit, so allowing some interest rate skew to be modelled and ensuring positive default 
intensities. Using the perturbation expansion techniques proposed by Horvath et al. (2017) we deduce 
an analytic expression for the relevant Green’s function operator which is asymptotically valid in the 
limit of low rates and intensities. We use this to produce approximate expressions for CDS prices. It 
is also shown how, converesly, these results can be used through a bootstrapping procedure to calibrate 
the model to a term structure of CDS rates. We consider also the pricing of defaultable Libor payments, 
which capability would be required for CDS with floating rate coupons, or else for interest rate swap 
extinguishers. 

1 Introduction 

We consider reduced form modelling of credit derivative prices where both the interest rate and the credit 
intensity (or spread) are considered to be stochastic and potentially correlated. We base our approach on 
that pioneered by Schonbucher (1999) who took both processes to be normal mean-reverting diffusions, in 
other words governed by the gaussian short rate model of Hull and White (1990). Solutions were found 
by constructing a two-dimensional tree. As was pointed out by Schonbucher (1999), it is a straightforward 
matter to extend his model to non-gaussian processes. A number of authors have followed this suggestion 
taking the credit process to be lognormal, governed by a Black and Karasinski (1991) short rate model 
which, although less tractable than a gaussian model, ensures that credit spreads stay positive (and thus 
that survival probabilities are decreasing functions of time). Jobst and Zenios (2001) sought to price portfolios 
of bonds, modelling the credit spread for securities in a given rating class in this way, coupled with a Hull- 
White interest rate model, but also allowing rating class migrations to take place. A similar approach with 
only rates and credit default risk was used by Cortina (2007) to provide analytic solutions for the prices of 
defaultable bonds in the assumed absence of correlation, and by Pan and Singleton (2007) who considered 
the joint distribution of credit spreads and default loss rates implied by CDS market data. 

We propose here a rates-credit hybrid model for CDS pricing. Interest rates are taken to be governed 
by the beta-blend short rate model introduced by Horvath et al. (2017): this is a one-parameter family of 
models where the choice of the parameter /3 allows the skewness to be fitted to the term structure implied by 
market prices of interest rate options. It provides an interpolation between the popular Hull-White (/? f 1) 
and Black-Karasinski (/3 = 0) models. The credit model we assume to be Black-Karasinski for the reasons 
already given. The two models we allow to be correlated, the impact of this correlation being our main 
interest in the present work. Horvath et al. (2017) has shown how these models can be solved analytically to 
provide Green’s functions which are asymptotically accurate in the limit of low rates. It is a straightforward 
matter to use these Green’s functions to calculate conditional values of future cash flows (in the interest rate 
case); likewise, in the case of credit modelling, we are able to obtain conditional survival probabilities and, 
by the same token, conditional default probability distributions. Our interest in the present work is to create 
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a bivariate Green’s function where both rates and credit are stochastic and possibly correlated. Typically 
it is only possible to price credit derivatives such as CDS in such circumstances numerically, but we show 
how, armed with an asymptotically valid representation of the requisite Green’s function, CDS prices can 
be inferred to good accuracy by analytic means from known credit spreads (or vice versa) in the presence of 
non-zero correlation. 

The layout of the paper is as follows. We set out our stochastic modelling assumptions and associated 
no-arbitrage conditions then derive the associated PDE in §2. We introduce in §3 an asymptotic scaling 
with the assumption that both interest rates and credit spreads are small, on which basis we are able to 
derive a Green’s function solution for our PDE as an asymptotic expansion. Calibration of the model to 
satisfy the specified no-arbitrage conditions is considered in §4. We use our Green’s function in §5 to derive 
the price of the protection leg of a CDS. It is shown in §6 how the (fixed) coupon leg is priced, taking into 
account the possible impact of accrued interest if this is required to be paid upon default. The reverse task 
of calibrating the model to CDS prices, if these rather than risky bond prices are considered to be known 
a priori, is considered in §7. We go 011 to show in §8 how defaultable Libor payments can be priced, which 
facility is required where a CDS (or credit linked note, or other credit-contingent instrument) pays floating 
rate coupons. We again take account of accrued interest payments at default. Finally, comparisons are made 
in §9 of results obtained using our formulae with more accurate finite difference results. They are seen to be 
highly favourable. 


2 Modelling Assumptions 

We model the interest short rate and the instantaneous credit default intensity as stochastic processes. These 
are taken, respectively, to be defined by a beta-blend model (see Horvath et ah, 2017) and the lognormal 
model of Black and Karasinski (1991). We shall find it convenient to work with auxiliary variables x t and 
y t satisfying the following Ornstein-Uhlenbeck processes: 

dxt = —a r Xtdt + a r (t) dW^, (1) 

dy t = - otxyt.dt + <7\{t) dWf, (2) 

where dW{ and dWf are correlated Brownian motions under the risk-neutral measure with 


corr (W t 1 ,W?)=p rX . 

These auxiliary variables are related to the interest short rate ry and the credit default intensity A*, 
respectively, by 


(1 - P)r t + Pr(t) = ( r(t ) + (1 - P)r*(t))£ , (3) 

X t = (X(t) + X*(t))£(y t ), (4) 

with p £ [0,1) assumed. Here f(t) is the instantaneous forward rate, X{t) the associated credit spread (see 
(8) below) and £(X t ) := exp (X t — \ [X]*) is a stochastic exponential with [X\ t the quadratic variation of a 
process X t } The required form of the configurable functions r*(t) and A *(t) is determined by calibration of 
the model to satisfy the 110 -arbitrage conditions set out below. We further assume that xq = yo = 0, with 
t = 0 the “as of” date for which the model is calibrated. As can be seen, the beta-blend model for the interest 
rate represents a hybrid between the Hull-White model (/3 f 1) and the Blaek-Karasinski model (p = 0). 
The credit intensity model is taken to be Black-Karasinski. 

1 We observe that, for /? £ (0,1), rt is not defined at times t for which r(t) = 0. This problem can be mitigated by replacing 
|r(t)|^ on the denominator in (3) with (|r(t)|+(l — /3)<5)' s for some small smoothing parameter S > 0. This will ensure a positive 
denominator as well as maintaining Hull-White and Black-Karasinski as the limiting case models when /3 —>- 1 and /? —► 0 
respectively. The analysis below goes through unaffected making this substitution throughout. 
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The no-arbitrage condition 

The formal no-arbitrage constraints which determine the functions r*(t) and A *(t) are as follows: 


E 



E 


g- /o(r»+A„)ds 


D(0,t), 


(5) 

( 6 ) 


under the martingale measure for 0 < t < T m , where T m is the longest maturity date for which the model is 
calibrated, 

D(t 1 ,t 2 )=e-ti? 7(s)ds (7) 

is the ti-forward price of the ^-maturity zero coupon bond and 

B(t 1 ,t 2 ) = e-£ ( 8 ) 


the corresponding risky bond price. We shall assume the bond prices can be ascertained at the initial time 
t = 0 from the market, whence we can view (7) and (8) as defining the forward rate r(t) and associated 
credit spread A (t), respectively. 


Derivation of governing PDE 

We consider the general problem of pricing a cash security with maturity T whose payoff depends on xt 
and jjT- We will also look in §5 at the case of a protection leg whose payoff depends on x T * and y r », where 
t* is a stopping time in (0,T]. We introduced the convenient shorthand notation that, for a process X t and 
deterministic function /(.,.), 

£x(f(X t ,t)) := £(f(X t ,t))\ Xt=x , 

in terms of which we can re-write (3) and (4) as r t = r(x t ,t) and A t = X(y t ,t), where 

r{x,t) := (W) + (1 - P)r*(t))£ & - pr{t)J , (9) 

\(y,t) := + (10) 

Writing the price of the security at time t £ [0,T] as fj = f(xt,yt,t), we can infer by application of 
the Feynman Kac theorem to (1) and (2) in the standard manner that the function f(x,y,t) satisfies the 
following backward diffusion equation: 

(J^ + £ - r(x, t) - A (y, t)j f{x, y, t) = 0, (11) 

where 

i := -“' i l - + l + *lm !$) (U) 

with in general frf = P(xT,yT ) for some payoff function P(.). 2 In the absence of closed form solutions to 
(11) and guided by the work of Hagan et al. (2005) and Horvath et al. (2017), we propose a perturbation 
expansion approach as follows. 

2 In the following we will not in fact consider payoffs dependent on y, although the methods employed could easily be extended 
to take account of such an eventuality. 
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3 Perturbation Expansion 

3.1 Definition of Scaled Variables 


For both short rate models we apply a "low rates” assumption. To this end we define, taking T m to be the 
longest time to maturity for which the model is calibrated, small parameters 


1 


€ r .— 


QL r Tm J 0 


r(t) dt 


£\ ■= 


a x T r 


m J 0 


and 0(1) functions 


r(t) := e. 


A (t) dt 

-i Ht) 
1-/3’ 




and 


Ht) = e A lA (t)> 

We further define a new scaled variable x t and associated scaled volatility a x (t) by 

x t := e~ ,3 Xte 0 ‘ rt , 

— u\„cx r t 


Likewise we define 


Tr (t) = e r p a r (t)e aT \ 


Vt ■= y t e axt , 
(T y {t) := <r x (t)e axt . 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


Here the exponential time scaling is to facilitate removal of the mean reverting drift terms in (12). We 
further define new functional forms /(•) and P(-) by: 

f(x,y,t) := B(p,t)fF\ Xt = Xtyt=y , 

P(x T ,yT ) = P(x T , yr), 


where xt and )/t are related to xt and yr by (17) and (18), respectively. In this notation, (1) and (2) can 
be re-expressed as 


dx t = cr x (t) dWf , 
dy t = cr v {t) dWf 

and (11) as 

(J^ + P ~ e r h(x, t) - e x g(y, t)j f(x , y,t) = 0 

where 

£| - 1 := l 

h(x, t) := h{x, t, t ) 
g(y,t) ■= g(y,t,t ) 


(19) 

( 20 ) 

( 21 ) 


( 22 ) 

(23) 

(24) 
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with 


h{x,t,t{) := (r(4) +r*(t 1 ))£ x {F 0 {t x )x t ) - r(t i), 4 > t, 

(\ - 

Fp{t) := U P) 


If(t )!/ 3 

g(y,t,h) := (A(t) + \*{t))£y (e~ axtl y t ) - A(4), 4 > t. 

We seek a Green’s function solution for (21) as a joint power series in e r and e\. 


(25) 

(26) 
(27) 


3.2 Derivation of Green’s Function 

We follow the methodology propounded by Hagan et al. (2005) and subsequently developed by Pagliarani 
et al. (2011) in observing that the operator 


r 

U(t',v)=ex p / (£{u) — e r h(x,u) — e\g(y,u))du 


(28) 


is a formal solution of (21) for t < v subject to U(v;v) = I. The Green’s function is then the integral kernel 
of li(t; v), viz. 

G(x,y,t;£,r),v) =U(t;v)(x,y,€,r)). (29) 

Note that the corresponding Green’s function for fj\ Xt - xyt - y will be given straightforwardly 3 by 


G*(x,y,t;£,ri,v) = B(t,v)G(x,y,t;£,r],v). 


(30) 


Again following Hagan et al. (2005) we consider first the limiting problem with e r — e\ = 0. The required 
solution can be written 

r v 

(31) 


r 

Go,o(t] v ) := exp / C(u)du. 


We will in all cases be interested in so-called free-boundary Green’s function solutions which tend to zero as 
x,y —> ±oo. The Green’s function solution subject to these conditions is well known. It is given by: 


<9 2 

Gq,o(x, y, t; r?, v) = ~x,y-y, R(t, v)), t < v 


(32) 


where N 2 (x,y; R(t, v)) is a bivariate Gaussian probability distribution function with mean 0 and covariance 
matrix 

4(M) 

I p (t,v) I y {t,v) 

with 


R(t, v) := 


(33) 


4(4,4) := / a~ x (u)du 
Jt ± 

ft2 


Pt2 

4(4,4):= / a'y(u)du. 

Jt-l 

ft 2 

4(4,4) ■= PrX / cr x (u)<7y(u) du. 
Jt 1 


(34) 

(35) 

(36) 


3 The required Green’s function must be of the form G*(.) = CB(0,t) 1 G(.). Correct normalisation of the new Green’s 
function at t = v then requires CB( 0, u)^ 1 = 1, from which we deduce (30). 
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Before deriving higher order terms in our proposed Green’s function expansion, we introduce some additional 
notation. First, following Horvath et al. (2017), we propose that the form of the configurable function r*(t) 
required to satisfy (5) is 

OO 

= (37) 

i= 0 

We further infer from the structure of (21) that the form of the expansion for A*(£) required to satisfy (6) is 

OO OO 

A*(<),yy<4Ah(t). (38) 

j= 0 i=0 

On this basis we can now expand h(x,t,t 1 ) as 

OO 

h{x , t, ti) = ^2 t, h) - f(ti) (39) 

i= 0 

where 

h 0 (x,t,ti) := (r(t i) + ro(ti))£ x (Fp(ti)x t ) 
hi(x,t,ti) := r*(ti)£ x (Fp(ti)x t ), i > 0. 

and similarly 

OO OO 

g{y,t,t i) ='^2'*Te' r e J x g iij {x,t,t 1 ) - \{t i) (40) 

j— 0 i—0 

where 

9 Qfi{y,tM) ■= (A(fi) + XQ 0 (ti))£y(e~ axtl y t ) 

:= X*j(ti)£ y (Fp(t 1 )y t ), i + j> 0. 

Note that the required quadratic variations are given by [x\t = /ic(0,f) and [y] t 

hi(x,t) := hi(x,t,t), i > 0, 

9i,j(y,t) ■= gi,j(y,t,t), i,j > o. 


= I y (0, t). We further define 

(41) 

(42) 


Asymptotic Expansion 

We propose that a formal asymptotic expansion of U(t;v ) is then possible as 

OO OO 

U(t;v) = 5ZX144 U i,j( t ’’ v ) ( 43 ) 

i—0 j—0 

where each of the Ui d {.) is 0(1). By a process of induction we infer that the following recurrence relation 
must hold: 

U id (t;v) = - ^ ^o,o(i;w) (ft(a:,ti)W i ^. liJ -(t;ti)l i> o + g{y,u)U id _i{u-,v)t j>0 )du, i,j > 0. (44) 

Substituting this expression into (43) and arid grouping like powers of e r and e\, we can write 

OO OO 

U(t\v) =J2J2e l r e{,U* d (t;v), 

i—0 j—0 
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where in particular we have, following Kato (1995), 

= U 0 ,o(t;v), 

Wi* 0 (t ; v ) = - j ^0,0 {t; u)(ho(x, u ) - r(u))U 0 ,o(u; v)du, 

^o,i (t]v) = - j Uo,o{t; «)(9o,o(9, u) - A (u))U 0 , 0 (u; v)du , 

W 2 * 0 (t; v) = - f Uofi(t; u ) ((/io(x, «) - f(u))Wi, 0 (u; v) + (/ii(ar, u) + 91,0(9, u))U 0 ,o(u] v)) du, 


Uo, 2 (t;v) = - / W 0 , 0 (t;u) ((5o,o(y, w) - A(u))W 0 ,i(u;w) + 90,1(9, u)U 0 fi (u\ v )) du 


v) = - J U 0 ,o (t; u) [ih 0 {x, u ) - f(u))U 0 ,i(u; v ) + (50,0(2/, «) - A(u))W 1]0 (u; v) 

+ 9i,o(5, u)U 0 ,o{u; v)j du. 

Taking Gjj(.) to be the integral kernel of U*j, we deduce: 

OO OO 

G*{x,y,t;£,T),v) = B(t,v)^2^2e l r e j x G i j(x,y,t;^,ri,v). (45) 

i=0 j= 0 

Following the approach taken by Pagliarani et al. (2011), we seek to obtain explicitly the form of the required 
integral kernels. At first order we obtain, after some manipulations: 

pv p 00 poo 

Gifi(x,y,t;£,r),v) = - / / / G 0 ,o(x,y,t;x 1 ,y 1 ,t 1 )(h 0 (x 1 ,t 1 ) - f(ti)) 

Jt j —OO J — OO 

Go,0(^1,51, h; £, 9, v) dx 1 dy 1 dt 1 

/ V 

(h 0 {x,t,ti)Ml tl - r(ti))Go,o {x,y,t;£,r),v) dt-i (46) 


and 


V P OO p OO 


G 0 ,i{x,y,t-^,r],v) = - 


Go,o(x, y, t\ xi,yi,ti)(go,o(yi,ti) - \{h)) 


t J —OO J —00 


G 0 ,o(xi, 51, h; £, 9, v) dx 1 dy 1 dt 1 


= ~ J (9o,o(9, - X(t 1 ))G 0 ,o{x,y,t;^,r],v)dt 1 , 

where we have defined the shift operators 

• A/ *ti,t 2 Go,o(2b9,^;£,9,«) : = G 0 , 0 (x,y,t;£- F / 3 (t 2 )I x (t 1 ,t 2 ),r],v) 
■M¥ 1 ,t 2 G 0 ,o(x,y,t;£,,ri,v) := G 0 ,o (x,y,t;^,rj - e~ axt 2 I v (t 1 ,t 2 ),v) 

Similarly at second order we obtain 

(h 0 (x,t,ti)Mt M -r(t{))(ho{x,t,t 2 )Mt lM -r(t 2 )) 

Go,o(x, y, t; £, 9 , v) dtidt 2 

PV pt 2 

+ (h 0 (x,t,t 2 ) (ho(x,t,ti)(exp(Fp(t 1 )Fp(t 2 )I x ( 0 ,ti))-l) 

Jt Jt 

9, £, 9,t>) dt x dt 2 

~ J (hi(x, t, h)M x t , tl + 91,0(9, t, ti)M y t , tl )G 0 ,o (x, y, t; £, 9, v) dt 1, 


(47) 

(48) 

(49) 


(50) 
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and 


Go,2(Z, y, t; t,y,v)= / (50,0(2/, t, ti)M v tM - \(ti)){g 0>0 (y, t, t 2 )M y tut2 - A (t 2 )) 

Jt Jt 

Go,o{x,y,t-,£,r),v) dtidt 2 

+ J* (50,0(2/, t,t 2 ) J (50,0(2/, Mi) (exp (e _aA(tl+t2) 4(0,ti)) - l) 

^ V tM M tut^ G ofi{x, 5, t\ £, 77, v) dhdt 2 

- J 9 o,i(yJJi)- M t,t 1 Gofi(x,y,t;£,ri,v)dt 1 ( 51 ) 

n t2 

(ho{x, t, ti)Mf tl - r(t 1 ))(g 0 , 0 (y,t,t 2 )M ! /^ t2 - \{t 2 )) 

Gofl(x,y,t;£,r),v) dt ± dt 2 

+ [ [ {9o,o{y,tM)M v tM ~\{ti)){ho(x,t,t 2 )M x tlM -f{t 2 )) 

Jt Jt 

Go,o(x,y,t-,t;,r],v) dt ± dt 2 

/ V r>t 2 

9o,o(x, t, t 2 ) J h 0 (x, t, ti) (exp ( e~ axt2 F l3 (t 1 )Ip(0 , tij) - l) 

Mt M M V t lM G 0fi {x, y, t; £, 77, v) dhdt 2 

rV r>t2 

+ ho{x,t,t 2 ) 50,0(5, Mi) (exp (e~ axtl Fp{t 2 )I p (S),ti)) - l) 

Jt Jt 

■Mt, tl Mt lyt2 G 0 fi(x, y, t; £, 77, v) dt x dt 2 

~ J 9i,o{yJJx)Ml tl G 0 ,o(x,y,t-,^r],v)dt 1 . ( 52 ) 

Substituting ( 46 )-( 52 ) into ( 45 ) gives the required Green’s function to second order accuracy. 


4 Calibration 

It remains to calibrate our model consistent with the no-arbitrage conditions (5) and (6). This is achieved 
by considering in the former case the consistent pricing of a risk-free cash flow, and in the latter case by 
going on to consider a risky cash flow, as we shall now show. 

Pricing of risk-free cash flow 

The calculation for a risk-free cash flow in our model is very similar to that performed by Horvath et al. 
(2017) and essentially corresponds to taking the distinguished limit as e\ —> 0 then e r —> 0. The same result 
is naturally obtained, namely that fj = X T (x,t ) where 

X T (x, t) ~ D(t, T) (l - e r F lfi (x, t) + efF 2 fi {x, t)) (53) 

with 0(e^.) errors, and our Green’s function gives rise to 

Fi,o(x,t) = [ (h 0 (x,t,ti) - r(ti))dt! 

hi(x, t, 1 1) dti 

+ J h 0 (x,t,t 2 ) J h 0 (x,t,ti)(exp(Fp(t 1 )Fp(t 2 )I x ( 0 ,ti))-l)dtidt 2 . 


2,o(x,t) = \F 2 0 {x,t) - Jt 


8 



Of interest, to us here is the conclusion that, setting x = y = t = 0 in ( 53 ), satisfying ( 5 ) above requires us 
to choose 


r* 0 (t) = 0 , ( 54 ) 

r i (i)=f(t) [ r(u){exp(Fp(u)Fp(t)I x ( 0 ,u))-l)du, ( 55 ) 

Jo 


Pricing of risky cash flow 

We continue by writing the price at time t, of a risky (zero recovery) cash flow at time T as / t T = Y T (x t , Vt, t), 
noting that, in this case, P(x, y) = 1 and /q = Y T ( 0 , 0 , 0 ) = B( 0 , T). We look to derive the general functional 
form of Y T (.) implied by our model, and in the process to determine the conditions on A *{t) necessary to 
satisfy ( 6 ) above. Applying our asymptotic Green’s function to this problem we conclude 

Y T (x,y,t) ~ B(t,T) (l - e r F lfi (x,t) - e x F 0 t i(y, t) + e^F 2 , 0 (a;,t) + e r e\F ltl (x, y, t) + e 2 x F oa (y,t)) ( 56 ) 

with 0 (e x + e\) error, where the Fi o(x,t) are as defined above and 


F o,i{y,t) 
F 0 ,2(y,t) 

Fi,i{x,y,t) 


■ = J t (50,0(2/,Mi) - A(ti))dti, 

: = h F o,i{y,t) - J 50 , 1 ( 2 /, t,t!)dti 

+ J 50,0(2/, t,t 2 ) J 50,0(2/, Mi) (exp (e -a>( * 1 + * 2 ) 4 ( 0 ,fi)) -l)dtidt 2 , 

:= F lt 0 (x,t)F 0 > 1 (y,t) - j 5 i,o( 2 /> t, h) dti 
+ ^ h 0 (x,t,t 2 ) J 50,0(5, Mi) (exp (e _a> * 1 F / 3(t 2 )/p(0,ti)) - l) dt ± dt 2 , 
+ / 50,0(5, M2) ^ h 0 (x,t,ti) (exp (e~ axt 2 F / s(t 1 )I p { 0 ,t 1 )) - l) dh dt 2 . 


Setting x = y = t = 0 , we find that the 110-arbitrage condition F T ( 0 , 0 , 0 ) = B( 0 ,T) is satisfied by the 
expression in ( 56 ) iff we choose 


= 0, 

Ao i(t) = A(t) J A (u) ^exp (^e~ ax< ' t+u ' > I y (0, u)^j — lj du. 

K 0 (t) = r(t) [ A (u) (exp (e~ axU F p (t)I p (0,u)) - l) du 

+ A (t) f r(u) (exp (eT axt Fp{u)I p {Q,u)) - l) du. 
Jo 

This completes the calibration of our model to second order. 


( 57 ) 

( 58 ) 


( 59 ) 


5 Pricing of Protection Leg 

We next extend our model to the pricing of a protection leg, viz. an instrument which pays 1 — R times the 
trade notional upon default of the named debt at some stopping time r*, zero otherwise, where R is the 
recovery level on the debt (assumed known in advance as is customary in CDS markets). Let us suppose 
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now that the process / f T = p(xt, yt , *) in our original unsealed notation represents the price per unit notional 
of an instrument offering protection until some future time T. This will satisfy: 

^ + £ ~ r{x, t) - X(x, tp(x u y t ,t) = -(1 - R)X(y, t). ( 60 ) 

subject to the final condition p(x,y,T) = 0. Switching to scaled notation and defining a function 

p(x,y,t) := e\(l — R) f t \ Xt — xyt — v , 

this can be seen to satisfy: 

(JJ t + C~ e r h(x, t) - e x g(y, t)j p{x, y, t) = -(A (*) + g(y, *)) ( 61 ) 

subject to the final condition p(x,y,T) = 0. As is evident, this is nothing other than a nonhomogeneous 
version of (21), hence can be solved (asymptotically) using the Green’s function G*(-) defined previously. 
Truncating the Green’s function in this case at first order and applying it in the standard manner, we deduce: 

p(x,y,t) ~ P 0 ,o(y,t) + e x Po,i{y,t) + e r P lt o(x,y,t) ( 62 ) 

with 0(e x + e\) error, where 

pT poo p oo 

£0,0(5,*)=/ / / B(t,t 1 )G 0 fi(x,y,t;^,r],t 1 )go,o{v,ti)d^dr]dt\ 

J t J —oo J —oo 
T 

£(*,*i)5o,o(2/>*>*i)<**i> (63) 



and 


pi pO o p OO 

£ 0 , 1 ( 2 /,*)=/ / / B(t,t 2 )Go,i(x,y,t;£,ri,t2)go,o{' r ht2) d£dr]dt 2 

J t J —00 J —OO 
pT poo poo 

+ / B(t,t 2 )G 0 fi(x,y,t;£,r),t 2 )go,i(ri,t 2 )d£dr]dt 2 

J t J — 00 J — 00 

= - ft £(*,*2)50,0(2/, *’*2)^ 5 o,o( 5 ,*,*i) (exp (e~ ax( - tl+i 2 ) I y (t,ti)^ -ljdhdt 2 
+ J t £(*,* 1 ) 50,1 (5, *, *i) d*i 

pT poo poo 

Pi,o{x> 5, *) = / / / £(*, * 2 )G f i,o(a;, 5, *; £, 5, * 2 ) 50 , 0 ( 5 , * 2 ) d£ dgdt 2 

J t J — 00 J — 00 


(64) 


ft J —OO J —00 
f-T poo poo 


't J —00 J —00 
cT 


£(*, * 2 )G 0 ,o(a:, 5, *; £, V, * 2 ) 51 , 0 ( 5 , * 2 ) <*£ di 7 d* 2 
*2 


/ -* f‘f'2 

£(*,* 2 ) 50 , 0 ( 5 , *,* 2 ) J ho(x, *, * 1 ) (exp (Fp(ti)e~ axt2 I p (t, * 1 )) — l) dt± dt 2 (65) 
+ ^ £(*,*1)51,0(5,t,*i)d*i- ( 66 ) 

Finally, noting that Po,i(0, 0) = 0 and combining the relevant contributions from £o,o(0,0) and Pi^O, 0,0), 
we conclude in terms of our original (unsealed) notation that the PV of the protection leg is given 
asymptotically by 


£Vt =0 ~ (1 - £) [ B(0,u) (X(u) (1 - <j>p(u)) + AAi(u)) dw 
do 


(67) 
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with 0(e\(e^. + e|)) error, where 


r(v) 

nt2 


ex P -p-TWTg - 1 I dv, 

kwr 


4 >p{u) : = 

:= P,A r e -(^+“A)(t2-«) fTr ( u)<TA ( u )) rf , 


AAi(tt) 

0i,i(«,f) 


0i,i(Wj u) dv, 


r(u)X(v) 

1-/3 


(1-/8 (O.v) 

ex P--1 

K«)r 


A(w)r(u) ^ (1 — /3)e ^7^(0, i>) n 

1 7 \ i n J- 


1 - , 


k(«)l^ 


( 68 ) 

(69) 

(70) 


(71) 


Note that AAi(t) here is none other than AJ 0 (t), expressed in unsealed notation. 

As expected, cancellation results in there being no contribution at O(e^). Further, in the absence of 
correlation, <f>p{-) = AAi(-) = 0 and the value of protection is as given under the assumption of deterministic 
rates. We further note that it is a straightforward matter to take the limit as /? f 1 in (68) and (71) to obtain 
the results for the case of a Hull-White interest rate model. In particular, we find 


0 p(m) 


0 —a r (u—v) 


Ir(0,v) dv 


AAr(u) -► J ( \{v)e- aAu ~ v) + A {u)e~ a ^ u - v ^ I R (0,v) 


dv 


6 Pricing of Fixed Coupon Leg 

For a payment period (t*_i, t»), i > 1, the payoff for a fixed coupon flow will be given by c6(ti-i,ti), with 
S(-) the year fraction under the relevant day count convention. The PV of this (defaultable) cash flow is 
given straightforwardly by discounting it with the relevant risky discount factor, namely 13(0, ti). If a default 
occurs at a time r* € (tj_i,tj), there may be a payment at time t* of the accrued interest amount, namely 

Accrual = C(5(t*_ i,r*), 


if such is contractually specified. 

We look to estimate the value of such payments under our model. Our task is simplified by realising that 
the value of accrual is essentially a forward-starting protection payment with a time-dependent payoff. The 
calculation is virtually identical to that in §5 above with the accrual amount incorporated into the r.h.s. of 
(61). We omit the details and simply state the result that, taking account of accrual, the value of the cash 
flow(s) associated with the ith payment period is given for U > 0 by 

rti 

PV t % ~ cB(0,ti)5(ti-i,U) + c / B(0,u)6(ti-i,u) (A (u) (1 - <t>p{y)) + AAi(u)) du (72) 

J 0Vij_i 

with 0(e\(ef, + e|)) error, and zero otherwise. Technically, since the coupon c should be considered to be 
0(e r ), the adjustment term here is, through the inclusion of <f>p(u), an order of approximation smaller than 
those already considered; however, we include it here for completeness. 

With the above results and those of the previous section, we are in a position to infer prices for credit 
default swaps (CDS) for a given assumed term structure of zero coupon bond prices, or alternatively to 
calibrate our model to a given term structure of CDS rates, as we shall now show. 
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7 Calibration to CDS Market 


In practice, our model would often be calibrated to a term structure of CDS prices rather than to zero coupon 
bond prices. We therefore seek to demonstrate how this could be done making use of the methodology 
developed in the preceding section. First of all, let us define the calibration CDS to have maturities given by 
Tj, j = 1,2,... ,n with corresponding par rates given by Sj. Let us further suppose that the function A(f) 
can be taken as piecewise constant between the Tj , given say by 

\(t) = Xj, t G (Tj-i,Tj]. 

We can write the PV of the coupon leg of the jth CDS with coupon amount c straightforwardly using the 
risky discount factors as 4 

Fj (c) = B(0,ti)5(ti-i,ti)l ti> o (73) 

i—1 

where to is the start date, the payment dates are t\,t2, ■ ■■, t„, t]\r j = Tj and S{.) specifies the relevant year 
fraction for the payment period. Denoting the PV of the respective protection legs by Pj, it is convenient 
to define the incremental protection values 


APj Pj Pj -1 
~ (i - R)\ :I 


B{ 0, u) (1 — <fip(u)) du 


(74) 


where, for a given estimate of A j, there will be a weak dependence of B(0,u) on that estimate, but no such 
dependence for Pj-i and we take Pq = 0. Successful calibration for the jth instrument requires 


A Pj = Fj(sj) — Fj-i(sj-i), 


with Fq = 0. These equations can easily be solved in turn iteratively by making an initial estimate, say 
A^ = si, or A^ = Xj-i for j > 1, computing Fj°\sj ) and A P^ based on the implied B(0,ii) for 
Tj_ i < u <Tj (the Fj_i(sj-i) term being already known from the previous calibration step) and using as 
the next estimate 



and so on until a converged value for A j emerges, for each value of j in turn. 


8 Defaultable Libor Flow 

Finally we give consideration to the pricing of defaultable Libor payments. This capability will be required 
in circumstances where a CDS is issued which pays floating rate coupons. This is particularly common for 
credit linked notes where, unlike with a CDS, an upfront notional is invested with protection sold thereon 
against a named debt issuer: typically Libor payments are made against the invested notional, plus a spread 
to take account of the credit risk. We consider a tenor-r credit-contingent, Libor flow paying at time T. This 
has zero value for t > T and can be treated as a fixed flow if t £ [T — t, T) and the Libor rate has already 
fixed. Hence we focus our attention on the case t < T — r. 

4 We assume here that there is no accrued coupon paid at default. The method to take into account the PV impact of 
contractual accrued coupon payments is described in §6 above. 
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Survival-contingent payment 

The payoff at time T of a tenor-r credit-contingent Libor flow is given in terms of our previously defined 
notation by 


Libor = X t (xt- t , T — r) 1 — 1 

= e r F 0 + D(T - r, T) -1 (e r F 1 } 0 (x T -r,T - r) + el(F 2 0 (x T - T , T - r) - F 2 fi (x T - T , T - r))) + O(e^), 


where 


F 0 := (1-/3) f r{t\) dt\. 

JT-t 

This final condition at time T requires specification of the unknown value of the stochastic variable Xt-t- 
We proceed by using our (first order) Green’s function first to value this payoff as of the fixing time t = T — r, 
at which Xt- t can be assumed known. The payoff can therefore be considered as a fixed (risky) cash flow 
and evaluated as such. Denoting the time-f value of the Libor flow by Fi,(xt,yt,t) and expanding as 

F l (x, y, t ) = e r H ly0 (x, t ) + e 2 r H 2 ,o{x, t) + e r exH 1 ^ 1 (x, y, t) + 0{e r {e 2 r + e 2 x )), (75) 

we obtain straightforwardly 


H h 0 (x, T-t) = B(T- t, T){F 0 + D(T - r, T)~ l F lfi (x, T - r)), 
H 2 ,o(x, T-t)= n { ^Z T ^] ( F lo(x, T-t)- F 2 , 0 (x, T - r)) 


- H lt 0 (x,y,T - t) [ (h 0 (x,T-T,ti)-r(ti))dti 
Jt-t 


IT- 

rT 


Hi,i(x,y,T - r) = -H lfi (x,y,T - r) / ( 50 , 0 ( 2 /, T - r, h) - X(h)) dh 


IT -T 


Taking this value as a payoff at time T — r and applying our leading order Green’s function for t < T — r, 
we obtain from Hi q(x,T — r): 


H\fi{x, t) = B(t, T) + D(T — t,T) 1 j (h 0 {x,t,h) - f(h)) dt-^j 


(76) 


13 



and consequently L/i,o(0,0) = B(0,T)Fq. In obtaining Hi ^(x,y,t), we observe that contributions will arise 
from applying Gq,i(-) to Hx,o(-) and from applying Go,o( - ) to H i j i(-)- We obtain: 


POO POO 

Hi,i(x,y,t) =-B(t,T- t) / G 0t i(x,y,t;£,r},T- t)H 1>0 (£,T - T)d£dr) 

J —GO J — OO 

pOO pOO 

+ B(t,T-r) / / G 0 fi(x,y,t;£,r],T -r)H 1A ((;,ri,T -T)d£dr] 

J —oo J —oo 
fT~r 

= —B(t, T) J {gofi{y,t,t 2 )~\{t 2 ))dt 2 

(f 0 + D(T - r, T)- 1 J ( h 0 {x , t, h) - f(ti)) dt^j 

f go,o(y,t,t 2 ) [ ho(x,t,ti) (exp (e- axt2 Fp(ti)I p (0,t. 2 )) — \) dt\dt 2 

JT-t 


D(T — t, T) J t 

- B{t,T) ^F 0 + D(T - T,r ) _1 j (ho(x,t,t 2 ) - r(t 2 )) dt^J J (ffo,o(y, *, h) - A(ii)) dt x 
B(t T ) 

"" ’ ' ho(x,t,t 2 ) / go,o{y,t,ti) 

Jt-t 

(exp (e~ axtl Fp(t 2 )I p (0, tx A t 2 )) - l) dtxdt 2 (77) 


D(T — r, T) 7 t _ t 


A similar expression can be derived for H 2 ^{x,t) by applying Gi.o(-) to iJ 10 (-) and Go o( - ) to H 20 (-) but, 
since in the absence of any e\ dependence it necessarily gives zero contribution at t = 0, we do not produce 
it explicitly here. Substituting (76) and (77) back into (75), setting x = y = t = 0 and reverting to unsealed 
notation, we see that our defaultable Libor flow has PV given by 


PV t=o ~B(0,T) 


1 - <pL 

D(T — t, T) 


-1 , 


with 0(e r (e^ + e^)) error, where 


rT - 


4*l ■= / x(v) 


f(„) j (i-»Wyy>A>) 


IT-r 


1 - 




J ' e - ax (v-u), u < v 

7 (u.v) := < ( , 

1 v ‘ ' j p—a. r (u—v) 


U > V. 


(78) 

(79) 

(80) 


Setting = 0, equivalent to zero correlation, recovers the price under the assumption of deterministic rates. 
As expected, negative correlation increases the value of the Libor flow (with risky discount factors B(0,t) 
assumed fixed) and vice versa. Also, taking the Hull-White limit as /3 f 1, we find 

v) / j(u, v)ljf (0, u A v) du dv (81) 

Jt-t 



We note that, provided A(-) remains bounded and a r ,a\ > 0, the value of <j>L in (81) remains bounded in 
the limit as T —» oo with r fixed, whence we conclude that the derived Hull-White formula is uniformly 
asymptotic in that limit. This will not be the case if either of the mean reversion rates are zero. Linder the 
assumption that f(-) also remains bounded, the same is seen to be true of (79). 


Adjustment for accrual 

For completeness we consider also the impact of any accrual payment which becomes due upon default. We 
focus on the accrual associated with a payment period (ti-x,U). The case ti-x > 0 can again be treated as 
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a fixed flow so we suppose w.l.o.g. that U-i > o* The payoff associated with default at a time t* £ 
of a credit-contingent Libor leg is given in terms of our previously defined notation by 

Accrual = A 1 (t*)(X u (x u _ 1 , U-i )^ 1 - 1 ) 

= e r A,(r*) (_Fg + 1 Fi 0 (x ti _ 1 , ti- 1)) + O(e^), 


where 


A / \ _ r, , I 

A»(' u) .— , ? wS [ij— 


( 82 ) 


d (ti—l , t j ) 

with S(-) the year fraction under the relevant day count convention. Writing the PV contribution associated 
with such a default as A H^\x,y,t) for t < ti and applying the same approach as in §5 above, we find the 
value as of time t £ is given by 


AH ( - l \x ti _ 11 y,t) = e r €x&HQ{x u _ x ,y,t) + 0 (e r e x (e r + e x )) 


( 83 ) 


where for t > ti -1 


fti poo poo 




B(t, v)G 0 fi (x, y, i; £, 77, 71)30,0(77, v)Ai(v) 


ft J — OO J — OO 


(F 0 + 1 F lfi (x,t l -i)) d^dr/dv, 


whence 


A Hl l J L (x,y,t i - 1 )= B(t i - 1 ,v)go,o{y,ti-i,v)A i (v)(F 0 + D(ti- 1 ,ti) 1 F lfi (x,U-i)) dv. ( 84 ) 

Jti -1 

Turning our attention to the case t < ti- 1, we take ( 84 ) as specifying the payoff at time fj_i; applying our 
leading order Green’s function we obtain straightforwardly: 

/ OO /»00 

/ B(t, ti-i)Go,o(x, y, t; £, 77, ti-i)Ai?$(£, 77, U-i) d£ dr] 

-00 J — OO 

rt 


= F 0 B(t,v)go l o(y,t i v)A i (v)dv + D(ti-i,ti) / B(t,v)g 0 ,o(y,t,v)Ai(v) 

Jti -1 Jti -1 

[ h 0 (x, t, u) (exp (e.- axV Fp(u)Ip{t, u A v)) — l) dudv, t < U- 1. ( 85 ) 


Setting x = y = t = 0 and reverting to original unsealed notation we obtain 


a pv: 


(i) 


t —0 


B( 0 ,v)X{v)Ai(v) 


1 + (j) A ( v) 


where 5 


u f{u) ( (1 — (d)'y(u, v)Ir( 0 , u Av) 


— 1 I dv. 


Mv)[= L,T^{ exp 


\r{u)\P 


— 1 du. 


( 86 ) 


( 87 ) 


Note that in the Hull-White limit (j f 1 we have 




7(77, v)Ir( 0 , u Av) du. 


'U-i 


( 88 ) 


5 As with the fixed coupon case, the accrual term is on order of approximation smaller than the main Libor contribution, 
so the <Pa{v) term here can be consistently ignored and is included only for completeness. However, if it is wished to include 
the impact of terms at 0(e\e (!) in the accrual term, arguably a term AAi(u) should also be included alongside A(v) in (86), 
capturing the impact of AJ 0 (-). 
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Figure 1: Impact of correlation on PV for a 5 y maturity CDS 


9 Comparison of Results 

A test is made of the accuracy of our CDS pricing formulae in the case of a Hull-White interest rate model. 
A 5y vanilla CDS was considered with a quarterly ATM coupon of 400 bp on a notional of 100M EUR. 
It was assumed that accrual was paid upon default. The zero rate was taken to rise from around 1.5% 
to around 3% over a five-year period. The credit spread was taken to be \{t) = 6.5% with an assumed 
recovery level of 40%. The (normal) interest rate volatility oy(f) was taken to be 0.5% with mean reversion 
rate a r = 0.25 and the (lognormal) credit volatility a\(t) was taken to be 60% with mean reversion rate 
a\ = 0.3. Calculations were carried out using our approximate formulae above and results compared with 
those from a finite difference solution of the relevant PDE. The impact of correlation on PV is illustrated 
for a range of correlations p r \ in Fig. 1. As can be seen, even up to a very large correlation of 80%, the 
differences nowhere exceed 0.002M EUR, viz. 0.2 bp of notional. In a similar calculation for a lOy CDS it 
was found that the maximum difference increased to 0.8 bp of notional, still of a negligible size. 

A similar test was done for an interest swap extinguisher. This is a trade where a Libor leg is exchanged 
with a fixed coupon leg, with all future cash flows being cancelled if a named issuer defaults. Typically, 
accrued coupon is paid up to the date of default. We consider the same fixed coupon leg as previously 
matched with a float leg paying Libor plus lOObp, which leaves the 5y swap quite close to the money. 
Calculations were again carried out using our approximate formulae above and results compared with those 
from a finite difference solution of the relevant PDE. The impact of correlation on PV is illustrated for 
a range of correlations p r \ in Fig. 2. As can be seen, the agreement is excellent. Notably the impact of 
correlation on the coupon flows is considerably less that on a protection leg. 
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